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Abstract 

Ion channels form pores across the lipid bilayer, 
selectively allowing inorganic ions to cross the 
membrane down their electrochemical gradient. 
While the study of ion desolvation free-energies 
have attracted much attention, the role of water in- 
side the pore is less clear. Here, molecular dynam- 
ics simulations of a reduced model of the KcsA 
selectivity filter indicate that the equilibrium posi- 
tion of Na+, but not of K+, is strongly influenced 
by confined water. The latter forms a stable com- 
plex with Na+, moving the equilibrium position 
of the ion to the plane of the backbone carbonyls. 
Almost at the centre of the binding site, the wa- 
ter molecule is trapped by favorable electrostatic 
interactions and backbone hydrogen-bonds. In the 
absence of confined water the equilibrium position 
of both Na+ and K+ is identical. Our observations 
strongly suggest a previously unnoticed active role 
of confined water in the selectivity mechanism of 
ion channels. 



1 Introduction 

Neurons enable us to think, act and remember.tn At 
the fundamental level, an important role is played 
by ion channels. Forming potassium-selective 
pores that span the cell membrane, potassium 
channels are the most widely distributed chan- 
nels in nature.!^ The breakthrough in the struc- 
ture determination came from the identification 
of the bacterial homolog from Streptomyces livi- 



r 



S2 
O23 



*To whom correspondence should be addressed 



t-; 



Figure 1: The KcsA channel. (Left) The whole 
protein channel. K+ ions are depicted in color 
(blue for the one in the S2 binding site). (Right) 
The reduced model of the S2 binding site, four 
di-glycines harmonically constrained to the crystal 
structure are used (see Methods for details). The 
experimental position of K+ at the center of the 
binding site, the extra and intra-cellular carbonyls 
are labeled as 5*2 , O12 and O23, respectively. 



dans (KcsA) .13 This channel is characterized by a 
tetrameric structure in which four identical protein 
subunits associate around a central ion conducting 
pore. At the extracellular side, the selectivity filter 
is formed by a highly conserved sequence of five 
residues (TVGYG).™ 

Originally, the high complementary of the selec- 
tivity filter to K+ was thought to be the reason for 
the selectivity. ^ But it is now clear that the inter- 
play between structure and dynamics of the selec- 
tivity filter is at the origin of the mechanism.!^ 
Recently, atomic models of the selectivity filter 
were used to elucidate the role of dynamics in the 
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process. '^^HEl in one of these models, the filter is 
reduced to the most selective binding site of KcsA, 
called S2.^ This is done by harmonically con- 
straining the four backbone segme nts d efining S2 
to the experimental conformation.li^^ Solvation 
free-energy calculations illustrated that the model 
is selective, allowing a statistical mechanics treat- 
ment of the limiting cases of rigid and very flexible 
binding sites. Moreover, calculations on larger 
models of the filter showed that backbone fluctua- 
tions are influenced by the presence of Na+ or K+ 
at position S2.^ These results support the view 
that several aspects of selectivity can be elucidated 
by analyzing ion binding to simplified models of 
the filter. 

Another important player in selectivity is wa- 
ter. Simulation results showed that confined wa- 
ter appears together with cations in the conduction 
pore.^'^^ Notwithstanding, it is not clear yet if 
water actively mediates selectivity or not. 

Here, the role of confined water is investigated 
by molecular dynamics simulations of a S2 bind- 
ing site model, providing evidence that the equi- 
librium position of Na+ within the binding site 
is displaced by the presence of a water molecule. 
Our calculations are in agreement with a recent 
crystallogra phic s tudy^and multi-ion free-energy 
calculations.l ^^l^^ l These concepts support the idea 
that KcsA can bind both Na+ and K+ with similar 
strength but different mechanism. 

2 Methods 

The S2 model. A reduced model of the S2 bind- 
ing site of the KcsA channel (PDB code: 1K4C, 
[1]) was built with four diglycine peptides as done 
in Ref. HD, The heavy atoms of the reduced 
model were constrained with an harmonic poten- 
tial of force constant k = 1000 kJ/mol/nm^ 2.4 
kcal/mol/A^). The coordinates were translated 
with the vector (—20,0,0), taking the axial coor- 
dinate Ap parallel to (1,0,0). A pdb file of the re- 
duced model is provided as Supplementary Infor- 
mation (SI). 

Molecular dynamics simulations. All calcula- 
tions were performed with the GROMA CS pr o- 
gramEQEU and the AMBER-03 force field.l^SESl a 
cubic box of initial length of 4 nm solvated with 



TIP4P-Ew water was used.f^ Simulations were in- 
tegrated with the Langevin equations (t = 0.2) 
at 300 K coupled with a Berendsen barostat 
(Tp = 1.0 ps).l22 Long range electrostatics was 
computed with PME^^ with a 1.0 nm cut-off for 
all non-bonded interactions. After 10 ns of equili- 
bration, each ion was simulated by a 100 ns long 
trajectory. For both Na+ and K+, the starting con- 
figuration was taken as the center of the binding 
site (52, see[T]). To check that there was no influ- 
ence on the starting position, 20 runs of 5 ns each 
were further performed (Figure SI in SI). Calcula- 
tions performed with a TIP3P water model are in 
agreement with the present analysis (see Figure S2 
in SI). 

Potential of mean force The potential of mean 
force for Na+ and K+ was computed with GRO- 
MACSl22l^ along the pore axis from Ap = —5.0 
(the bulk) to A^^ = 0.1 (position S2). Umbrella 
sampling calculations were spaced by 0.1 A along 
the axial coordinate Ap with the ion restrained in 
the normal plane with an harmonic potential (ky = 
/:^=1000 kJ/mol/nm^). An additional harmonic po- 
tential of force constant = 10000 and 20000 
kJ/mol/nm^ was applied in the —5.0<x< —2.1 
and —2.0 < x < —0.1 range, respectively. After 
1 ns of equilibration, trajectories were run for 1 ns. 
The weighted histogram method was used to re- 
construct the potential of mean force. 

Ion interaction energy in vacuo. To cal- 
culate the ion interaction energy Eion, the 
cation was harmonically restrained {{kx^ky^k^) — 
(50000, 1000, 1000) kJ/mol/nm^) at 0.25 A spaced 
positions along the axial coordinate Ap from -5.0 
to 0.0 A. After 1 ns of equilibration, each run was 
performed for 10 ns at 300 K and constant volume. 

3 Results 

The KcsA channel is shown in [T] In our simu- 
lation study, a reduced model of the protein se- 
lectivity filter was used. The S2 binding site was 
modeled by four peptides constrained to the ex- 
perime ntal st ructure (right panel, see Methods for 
details) .fi^l^ Within this model, the pore axis cen- 
tered along the channel is labeled as Ap, while the 
origin of the axis is taken as the experimental posi- 
tion of the K+ ion. This position is conventionally 
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Figure 2: Probability density functions along 
the pore axial coordinate Ap. (Top) K+ (black 
countoured blue area); (Bottom) Na+ (black coun- 
toured orange area). Position of the closest water 
molecule and ofOu and 023 carbonyls are showed 
as blue areas and dashed lines, respectively. 



called 52. Similarly, the positions of the carbonyl 
oxygens defining the binding site at the extra and 
intra-cellular sides of the model are denoted as Ou 
and O23, respectively ([!]). 

The behavior of and Na+ inside the S2 bind- 
ing site was studied by molecular dynamics simu- 
lations in explicit water (see Methods for details). 
In both cases, the starting position of the ion was 
S2. |2] shows the probability distribution of the ion 
position on the Ap axis. As expected, was 
found at position 5*2 (black countoured blue area) 
between Ou and O23 (dashed lines), coordinat- 
ing with the eight carbonyl oxygens of the binding 

This is not the case for Na+. After few ns of 
simulation (see Figure SI in SI), the ion hopped 




Figure 3: Stable configurations inside the KcsA 

filter. (Top) Na+ is in-plane with the O23 car- 
bonyls. The presence of a confined water molecule 
inside the filter stabilizes this position. As shown 
by the 90 degrees rotated structure, this water 
molecule makes alternate hydrogen bonds (red 
thick lines) with all four Ou carbonyls. (Bot- 
tom) K+ is stably bound into the conventional 5*2 
binding site. Outside the filter, one and two wa- 
ter molecules are present for Na+ and K+, respec- 
tively (ghost waters). 

from 5*2 to 023, assuming a configuration perfectly 
in plane with the 4 carbonyl oxygens (black coun- 
toured orange area in bottom panel of [2]). This po- 
sition is very stable, representing the 93.6% of the 
total simulation time, with a high barrier to hop 
back to 5*2 (roughly four transitions in 100 ns). 
Analysis of the closest water molecules to the ion 
provided a mechanism for the configuration shift. 
Contrary to the case of K+ where the solvent is 
at the outside of the pore, one water molecule en- 
ters the channel effectively shifting the position of 
Na+ to 023- This is shown by the probability dis- 
tribution of the closest water to the ion represented 
as a light blue area in the figure. A structural rep- 
resentation of the confined water is illustrated in 
|3| Per se, there is no energetic preference in shift- 
ing Na+ to the 023 position as shown by the aver- 
age interaction energy AEion between the ion and 
the pore in vacuo (|4]). In the absence of water, 5*2 
is the most stable position for both K+ and Na+, 
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Figure 4: Average interaction energy in vacuo 

between the ion and the binding site. Data for Na+ 
or K+ is shown in orange and blue, respectively. 
The position of the O23 oxygen carbonyls is shown 
as a vertical dashed line. 



strongly indicating that the binding site shift is due 
to the presence of a confined water. 

The confined water at position 5*2 is stabilized 
by multiple strong contacts ([3]). The water oxygen 
makes electrostatic interactions with Na+, while 
hydrogen bonds with two of the four carbonyls at 
position O12 are formed ([3} top right panel). These 
hydrogen bonds are also entropically stabilized, 
being the confined water capable of binding to all 
four carbonyls of O23 by rotating itself around the 
Ap axis (bond lifetime around 2 ps, see Figure S3 
in SI). This molecule is extremely stable, never ex- 
changing with the bulk. At the outside of the chan- 
nel another water molecule was found at the other 
side of Na+, forming favorable electrostatic inter- 
actions with the ion (ghost water in the top panel 
of |3]). The position of the confined water correlates 
with the fluctuations of Na+ around 023. indicat- 
ing that the ion- water configuration behaves as a 
complex ([5]). In the rare occasions when Na+ hops 
back to S2, the water molecule is expelled and the 
complex broken (Ajf^ ^ and A^ ^ S.sM, 

On the other hand, for the case of water 
was only found at the outside of the pore (ghost 
molecules at the bottom of [3]). Those waters are 
stabilized by the formation of hydrogen bonds 
with the carbonyls of the binding site. But the in- 
teraction with the ion is much weaker in this case, 
being the water oxygens facing the bulk. 

To complement the simulation study, free- 
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Figure 5: 2D Probability density function of the 
ion (Ap^) and confined water (Ap positions for 
the case of Na+. The contour lines are drawn for 
the probability values: 3 • 10"^ 10"^, lO-\ 1.0, 
4.0 and 8.0. 
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Figure 6: Potential of mean force for K+ and Na+ 

along the axial coordinate Ap (thick lines). The 
probability Py^ to have a water molecule inside the 
filter is represented as thin lines. Data for Na+ or 
K+ is shown in orange and blue, respectively. 
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energy calculations were performed. In [6] the po- 
tential of mean force (PMF, see Methods) along 
the axial coordinate Ap is shown (thick lines). 
The most stable configurations for Na+ and K+ 
were respectively found at positions O23 and 5*2, 
confirming our analysis. For both cases, posi- 
tion 023 is coupled with the presence of a water 
molecule inside the filter {P^ ^ 1 around 023. thin 
lines). Though energetically unfavorable, K+ at 
this position is in complex with a confined wa- 
ter as observed for Na+. Interestingly, the free- 
energy difference between the outside of the pore 
(Ap = — 5) and the most stable binding position 
is remarkably similar for the two ions, indicating 
no strong preference towards K+ binding. These 
results strongly support the idea that the filter has 
the ability to bind both ions with different mecha- 
nisms but similar strength as already suggested in 
Ref.18 1^ 



librium configurations. 

Our calculations on a minimalistic model of the 
filter provided a mechanism for the position shift. 
The presence of a confined water molecule in com- 
plex with Na+ effectively modifies the equilibrium 
configuration. The confined water is stabilized 
by favorable electrostatic interactions with the ion 
as well as multiple hydrogen-bonds with the Ou 
backbone carbonyls. The binding position shift 
disappears in the absence of confined water, mak- 
ing the latter an essential ingredient for the pref- 
erential position of Na+. In-plane binding always 
implies the presence of confined water even in the 
unfavorable case of K+. 

In the past, water was mostly considered in terms 
of solvation free-energies and its screening effects 
without much attention to the molecular mecha- 
nism. Our results reinforce the idea that biological 
water is an active player at the molecular level. ^ 



4 Discussion 

It has been known, for more than a decade now, 
that the selectivity filter of KcsA presents sev- 
eral binding sites. Each of them forms a cage of 
optimal coordination for K+ by means of eight 
backbone carbonyls. Most of the calculations on 
KcsA selectivity were based on the relative stabil- 
ity of Na+ over K+ inside the pore with respect 
to the bulk.l^^^ That is, the solvation free-energy 
to move an ion from bulk water to a specific bind- 
ing site inside the selectivity filter. These works 
suggested the S2 binding site as the most selective 
portion of the filter, being the free-energy differ- 
ence in the 5*2 position much more favorable for 
K+ compared to Na+. But this approach is not 
without problems if S2 is not a stable configura- 
tion for Na+. In fact, Na+ and K+ might be char- 
acterized by distinct equilibrium positions inside 
the pore. X-ray crystal logra phy*^ and multi-ion 
free-energy calculations^^ supported this idea, 
showing that Na+ preferentially adopts a configu- 
ration in-plane with backbone carbonyls. The ob- 
servation of distinct binding positions has several 
consequences to our understanding of ion selectiv- 
ity. Previous free-energy calculations using posi- 
tion restrains to the center of the binding sites need 
to be extended taking into account the correct equi- 
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